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TWO  DIMENSIONAL  FLUID  SIMULATION  OF  RELATIVISTIC 
ELECTRON  BEAM-HIGH  Z TARGET  INTERACTIONS 


I.  INTRODUCTION 

Calculations  of  pinched  flows  in  relativistic  electron  beam 
diodes  have  so  far  relied  heavily  on  particle  simulation  codes.1-3  How- 
ever, another  possible  approach  consists  of  treating  the  electron  beam 
as  a fluid  along  with  the  anode  plasma.  Both  approaches  might  be  com- 
plementary in  the  understanding  of  mechanisms  taking  place  inside  the 
diodes . 

In  particular.  Refs.  2 and  3 have  been  highly  successful  in  ex- 
plaining the  pinching  of  the  electron  beam  by  bootstrapping  its  way 
toward  the  axis  of  the  diode,  showing  how  anode  ions  and  plasma  aid  in 
the  pinch  formation.  However,  these  works  make  specific,  simple  assump- 
tions concerning  the  electromagnetic  fields  and  the  motion  of  the  anode 
plasma. 

It  is  clear  that  a fully  self-consistent  two-dimensional  model  of 

a relativistic  electron  beam-plasma  interaction  is  desired.  In  this 

present  work,  we  report  on  the  numerical  study  of  such  an  interaction  by 

means  of  a two-dimensional  fluid  model.  In  this  first  attempt,  we  have 

made  simplifying  assumptions  (convenient  for  a fluid  model)  and  in 

particular,  it  is  assumed  that  the  plasma  is  low  enough  in  density  to 

neglect  beam-plasma  collisions  so  that  the  electrons  can  be  treated  as 
Note:  Manuscript  submitted  April  12,  1977. 
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cold.  Calculations  start  when  the  self-field  of  the  beam  starts  to 
diffuse  into  the  plasma.  It  is  desired  to  determine  how  charge-neutrali- 
zation provided  by  the  plasma  causes  the  incident  beam  to  pinch  in  the 
diffused  magnetic  field. 

In  Sec.  II,  a general  model  for  the  relativistic  e-beam  plasma 
interaction  is  presented  as  well  as  initial  and  boundary  conditions. 

The  general  model  includes  collisions  of  the  beam  with  the  plasma.  Sim- 
plified equations  solved  in  the  remainder  of  this  study  are  then  pre- 
sented and  their  limits  of  application  discussed.  In  Sec.  Ill,  the 
computational  results  obtained  are  displayed.  In  Sec.  IV,  conclusions 
are  summarized  and  further  steps  are  suggested  in  order  to  solve  the  full 
set  of  general  equations. 

II.  MODEL  EQUATIONS 

Relativistic,  collisional  electron  beam-plasma  equations  have  been 
derived  and  reported  by  Mosher.4  They  are  presented  here  again  for  the 
sake  of  completeness.  If  n denotes  the  beam  electron  density,  ju  its 
fluid-flow  velocity,  ( v ) the  value  of  y (the  relativistic  factor) 
averaged  over  the  beam  electron  distribution  function,  Oi and  Oh  the 
effective  values  of  Z A and  Z2  ir.  A for  interaction  with  plasma 
electrons  and  ions  respectively,  then  the  beam  fluid  equations  can  be 


written: 
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where  n^  is  the  plasma  ion  density,  v the  plasma  fluid  velocity,  p the 
plasma  pressure  (=  (1  + Z)n^kT)  , £ the  total  plasma  energy  density, 

and  j the  plasma  current.  The  following  Maxwell  equations  close  the 
above  system  when  quasineutrality  is  assumed: 


_ BB 

7 x E , 
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7 x B = 4tt(  1 - enu)  . 
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The  system  of  Eqs.  (l)  - (3)  describes  the  full  interaction  of  a 
relativistic  electron  beam  with  a plasma  once  the  plasma-electron  trans- 
port equations  are  included.  Let  us  note  at  this  point  that  the  time- 
scales  for  the  electron  and  ion  dynamics  are  very  different  (of  the  order 
of  3 orders  of  magnitude)  and  constitute  one  of  the  major  problems  to 
the  simultaneous  solutions  of  these  equations.  Dealing  with  both  scales 
in  the  same  computer  run  would  be  costly  and  wasteful  unless  a quasi- 
static approach  to  solving  the  electron  equations  is  attempted  in  order 
to  eliminate  the  short  time- step  requirement  imposed  by  the  Courant  con- 
dition. This  is  reasonable  since  relativistic  electron  transit  times 
are  much  smaller  than  the  time  over  which  macroscopic  quantities  vary. 
Such  a formulation  is  in  development.  Equations  (1)  - (8)  can  be 


r 


t . .1  » n.i  , 


' v.. 


v-^arasr  se*r 


simplified  depending  on  the  region  of  the  beam  plasma  interaction  under 

investigation.  In  fact,  the  electron  beam-anode  plasma  system  can  be 
separated  into  three  regions.  In  the  first  region,  the  beam  flows 
through  a vacuum  containing  field- emitted  ions  from  the  anode3  and  is 
governed  by  the  action  of  space-charge  electric  and  magnetic  fields  and 
the  applied  axial  accelerating  field.  In  this  region,  Eqs.  (l)  - (j5) 
apply  with  the  coupling  term  set  equal  to  0 and  E determined  from 
Poisson's  equation.  In  Region  II,  the  beam  interacts  collectively  with 
a low  density  plasma  with  diffusion  of  the  magnetic  field  determined  by 
the  plasma  conductivity  in  that  region.  In  Region  III,  the  beam  inter- 
acts collisionally  with  the  high-density  plasma  produced  near  the  anode 
surface.  The  full  system  of  Eqs.  ( 1)  - (8)  apply  in  that  region  where 
collisions  dominate.  However,  the  action  of  dynamic  friction  which 
broadens  the  beam  energy  distribution  is  difficult  to  determine  without 
recourse  to  kinetic  theory.0*7 

In  this  work,  we  limit  consideration  to  Region  II  where  both 
quasineutrality  and  no  collisions  can  be  assumed.  In  that  region  where 
collisions  are  not  important  but  where  the  field  diffuses,  the  full 
system  of  equations  can  be  reduced  to  the  following: 
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7 x B = W ( j - enu) 


In  these  equations,  it  has  been  assumed  that  the  ion  dynamics  is  frozen 


and  only  the  beam  equations  and  Maxwell  equations  are  solved  for.  The 


plasma  conditions  enter  these  equations  through  an  assumed  constant 


scalar  conductivity 


j = ctE  . 


Furthermore,  it  is  assumed  that  the  electron  beam  is  cold  (i.e.  the 


pressure  term  is  equal  to  0)  since  the  beam-plasma  collisions  are  negli- 


gible. As  a result,  the  energy  of  the  beam  shows  only  in  kinetic  form 


and  remains  unchanged  fi.e.  changes  in  electric  potential  inside  the 


plasma  « 
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so  that  the  beam  energy  equation  is  not  required  any 


more.  Equations  ( 11)  - (12)  can  be  shown  to  constitute  the  diffusion 


equation  for  the  magnetic  field.  In  Eq.  (12),  = - enu  is  the  primary 


current  and  j = crE  is  the  plasma  current.  Note  that  when  = - j 


the  total  current  is  0 and  we  have 


7 X B 3 0 . 


The  initial  and  boundary  conditions  are  the  following.  At  t = 0 


and  z = 0 (which  delimits  the  beginning  of  Region  II),  the  beam  front 


enters  the  low  density  anode  plasma.  In  that  plane,  the  magnetic  field 
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The  initial  electron  beam  density  distribution  in  the  radial  direction 
is  chosen  as  a square  profile  whose  edge  has  been  smoothed  slightly  for 

computational  purposes.  The  vacuum  density  outside  the  beam  has  been 
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chosen  to  be  10  times  the  initial  beam  density  for  the  present  calcula- 
tions, but  it  can  be  changed  easily.  The  radial  velocity  u^  has  been 
taken  equal  to  0 everywhere  at  t = 0.  The  radial  electric  field  is 
assumed  to  be  0 at  z = L where  L delimits  the  end  of  Region  II.  This 
condition  results  from  the  assumption  that  the  anode  is  an  equipotential 
plane. 


Another  characteristic  of  the  model  is  that  a,  the  only  quantity 
characterizing  the  plasma  in  the  equations  solved,  is  assumed  to  be  con- 
stant, corresponding  to  an  isothermal  plasma  in  Region  II.  For  practical 
reasons,  c was  chosen  so  that  complete  diffusion  of  the  magnetic  field 
occurs  after  a few  hundreds  of  time  steps.  This  requirement  corre- 
sponded to  a very  low  temperature  of  the  plasma  ( 10  eV  and  below).  In 
reality,  <7  is  determined  from  the  plasma  temperature  achieved  by  electron- 
beam  heating  in  Eq.  (6). 

Numerically,  the  grid  has  been  chosen  inhomogeneous  in  the  r- 
direction  in  order  to  increase  the  numerical  resolution  near  the  axis. 

The  spatial  step  goes  from  0.1  cm  near  the  axis  to  3.6  cm  in  the  vacuum. 
Equations  (9)  and  ( 10)  are  solved  explicitly  using  FCTe  whereas  the 
magnetic  field  diffusion  equation  is  solved  implicitly.  The  time-step 
is  controlled  by  the  Courant  condition  based  on  the  electron  fluid 
motion. 


III.  RESULTS 

The  case  presented  below  corresponds  to  the  following  input  param- 
eters : 


Beam  radius 
Beam  current 
Beam  axial  velocity 
Length  of  Region  II 


5 cm 
1 MA 

v 

S2  - — - 0.9^3  (y  = 3) 

c 

0.5  cm 
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Figures  1 and  2 show  density  contours  at  t = t^^/2  and  t = 2 t 

respectively  where 


diff 


_ b-rjcL2 

Tdif  f 2 ' 

c 

The  initial  density  was  set  equal  to  2.8  x lO1^  and  remains  so  in 
the  plane  z = 0.  Figure  1 shows  pinch  formation  occurring  obliquely 
with  respect  to  the  diode  axis.  The  maximum  density  at  the  time  of 
Fig.  1 is  2 x 1012  , almost  an  order  of  magnitude  above  the  initial 
density.  Figure  2 shows  strong  pinching  on  the  axis  near  z = L.  The 
density  is  then  almost  200  times  the  initial  density. 

Figures  3 and  k-  show  e-beam  fluid  streak  lines  at  the  same  times 
as  Figs.  1 and  2.  The  pinching  is  seen  to  form  closer  to  the  end  of 
Region  II  as  the  radius  decreases.  As  time  develops,  this  remains  true 
but  the  pinching  recedes  toward  the  cathode.  Outside  the  pinching 
region,  we  note  areas  where  deconfining  forces  dominate  temporarily. 

The  magnetic  field  at  r = 3.21  cm  is  shown  in  Fig.  5.  It  is 
clear  from  this  graph  that  the  field  excluded  initially  progressively 
fills  the  plasma  in  Region  II.  By  t = 2 the  field  has  diffused 

entirely  into  the  anode  plasma. 

Figures  6 and  7 show  radial  profiles  of  the  magnetic  field  at 

t ■ t,.j././2  and  t =*  2 At  z = 0,  the  self-field  of  the  beam  is 

diff  diff 

shown,  increasing  linearly  with  r up  to  the  beam  boundary  and  then  de- 
caying as  Vr  in  the  vacuum.  Departures  from  the  linear  rise  indicates 
a higher  current  density  near  the  axis,  corresponding  to  pinching.  This 
pinching  is  shown  very  clearly  at  z = L in  Fig.  7.  The  beam  axial 
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velocity  is  down  by  more  than  a factor  of  10  in  this  area.  One  of  the 
results  found  with  our  model  is  that  no  steady  state  has  been  reached 
after  2 and  the  flow  remains  unsteady  in  all  parts  of  the  Region  II 

of  the  diode  under  consideration. 

IV.  CONCLUSIONS 

In  this  work,  we  have  described  the  electromagnetic  interaction 
of  a relativistic  e-beam  with  an  anode  plasma.  We  have  shown  the  solu- 
tion in  the  low-density  plasma  region  wnere  fields  diffuse  but  where 
collisions  are  unimportant.  We  have  found  that  the  beam  pinches  in  such 
a configuration  but  that  the  solution  does  not  reach  a steady-state  for 
the  time- scales  involved  in  this  study.  This  is  in  part  due  to  the  fact 
that  the  plasma  does  not  respond  to  beam-heating  so  that  magnetic  diffu- 
sion continues  to  larger  values  of  z.  If  the  plasma  was  allowed  to  heat, 
high  temperature  in  the  pinch  would  stop  the  diffusion.  In  any  case, 
the  pinching  observed  in  the  regime  under  investigation  is  significant 
because  it  shows  that  a low  density  plasma  in  front  of  the  anode  can 
cause  its  formation. 

The  discussed  solution  constitutes  a first  step  in  the  solution 
of  a full  beam-plasma  dynamic  interaction.  The  relativistic  e-beam 
fluid  equations  have  been  tested  as  well  as  the  magnetic  field  diffusion 
equation.  Addition  of  a pressure  term  to  the  beam  equations  and  colli- 
sional  coupling  to  a plasma  system  of  equations  remain  to  be  done.  In 
view  of  the  length  of  the  calculations  required  for  the  e-beam  dynamics 
when  performed  on  an  electron  time-scale,  this  study  strongly  suggests 
that  the  beam  equations  be  solved  implicitly  if  coupling  to  a much  slower 


dynamics  is  planned  as  a next  step.  Nevertheless,  this  study  has 
provided  interesting  results  of  its  own  and  pointed  out  causes  of 
eventual  problems  in  an  implicit  treatment. 
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Fig.  5 — Magnetic  field  intensity  at  r = 3.21  cm  as  a function  of  time 
Note  initial  conditions  and  final  condition  (t  = 2 


Fig.  7 — Magnetic  field  radial  profiles  at  t * 2 rdiff  at  various  z. 
Note  that  at  z * 0,  the  B field  radial  profile  does  not  change. 
Also,  large  values  of  field  near  axis  indicate  pinching. 
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